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The charge carrier drift mobility in disordered semiconductors is commonly determined from a single transit 
time graphically extracted from time-of-flight (ToF) photocurrent transients. However, the term transit time 
is ambiguously defined and fails to deliver a mobility in terms of a statistical average. Here, we introduce an 
advanced computational procedure to evaluate ToF transients, which allows to extract the whole distribution of 
transit times and mobilities from the photocurrent transient, instead of a single value. This method, extending 
the work of Scott etal. (Phys. Rev. B 46, 8603), is applicable to disordered systems with a Gaussian density of 
states (DOS) and its accuracy is validated using one-dimensional Monte Carlo simulations. We demonstrate the 
superiority of this new approach by comparing it to the common geometrical analysis of hole ToF transients 
measured on poly(3-hexyl thiophene-2,5-diyl) (P3HT). The extracted distributions provide access to a very 
detailed and accurate analysis of the charge carrier transport. For instance, not only the mobility given by the 
mean transit time, but also the mean mobility — which are generally not identical — can be calculated. Whereas 
the latter determines the macroscopic photocurrent, the former is relevant for an accurate determination of the 
energetic disorder parameter a within the Gaussian disorder model (GDM). a derived by using the common 
geometrical method is, as we show, underestimated instead. 



I. INTRODUCTION 

Disordered or amorphous semiconductors provide the po- 
tential of cheap future electronics with a wide variance of 
physical properties, e.g., absorption spectrum. This flexibility 
allows to design them for specific areas of application, rang- 
ing from transistors and light emitting diodes to sunlight har- 
vesting and converting solar cells. A common key feature de- 
termining the electrical device performance in these materials 
is the hopping type charge carrier transport. Therefore, a lot 
of effort has been put into studies during the last decades to 
understand the detailed mechanism, experimentally and theo- 
retically. 1 " 6 

In classical descriptions the charge transport is based on the 
electrostatic force induced drift and diffusion due to a spatial 
charge carrier density gradient. The related parameters are the 
charge carrier mobility /j and charge carrier diffusivity or dif- 
fusion constant D, respectively, which are connected by the 
classical Einstein 7 -Smoluchowski 8 relation D//j = kgT je. 
Here, kg is the Boltzmann coefficient, e the elementary charge 
and T the temperature. This concept successfully describes 
semiconductor devices coupling continuity and Poisson equa- 
tions 9 ' 10 in combination with physical models concerning, 
e.g., dissociation 11 ' 12 or recombination. 13 ' 14 

However, in the case of disordered systems it is known that 
the concept of mobility can be ill-defined 15 or, in other words, 
the transport is determined by a broad distribution of mobil- 
ities. This fact is due to the random nature of the charge 
carrier transport in an energetically and spatially disordered 
environment. It leads to energetic relaxation of the charge 
carriers with time, 15 " 18 and to an inhomogeneous distribution 
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of transport properties throughout the thin-film device. A 
detailed experimental investigation of the relation between the 
macroscopic behaviour and the microscopic transport statis- 
tics, however, requires the extraction of the transport param- 
eter distributions. Furthermore, new conclusions in terms of 
the energy distribution of the DOS may be drawn from the 
shape of these distributions. Such an approach should signif- 
icantly improve the ability to check the consistency between 
experimental data and various theoretical models. In turn, we 
consider it as a fundamental step towards a more general func- 
tional description of the charge carrier transport in disordered 
systems in terms of drift and diffusion. In some cases, such 
descriptions are already available, 15 ' 16 ' 18 ' 21 ' 23 " 26 but not in the 
general form presented here. 

First attempts to experimentally extract mobility distribu- 
tions from time-of-flight (ToF) 27 " 29 measurements or from 
the turn-on dynamics of polymer photo-cells where published 
by Scott, etal? and Rappaport etal. 20 ^ 1 , respectively. Un- 
fortunately, the method by Scott etal. was seldom used since 
its publication in 1992. 32 The same holds for the method 
published by Rappaport etal. Moreover the full potential util- 
ising the statistics of the extracted parameter distributions 
was hardly ever discussed. In the case of Scott's method 
this is probably due to the fact that the method as such is 
not suitable for systems exhibiting a Gaussian density of 
states (DOS), found in many organic semiconductors, e.g., the 
widely used conjugated polymer poly(3-hexyl thiophene-2,5- 
diyl)(P3HT). 

In this paper we extend the method by Scott et al. with re- 
spect to a Gaussian DOS and reveal the advantages of this ap- 
proach over the conventional geometrical technique to extract 
mobilities. The geometrical method determines the charge 
carrier mobility from a particular transit time, which is graph- 
ically extracted from ToF transients. This transit time is 
defined by the intersection of the fits to the pre- and the 
post-transit sections of the photo current transient, respec- 
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lively • . Henceforth we refer to the methods compared in 
this work as the Scott method and the geometric method, re- 
spectively. To verify the accuracy of the our new approach 
and the modifications to Scott's method we utilise a one- 
dimensional Monte Carlo simulation. Generally, the justifi- 
cation of an alternative, new approach is best demonstrated in 
comparison to a well-established standard method. Hence, the 
comparison is done here by evaluating the hole ToF transients 
measured on P3HT with the Scott and the geometric method, 
respectively. Subsequently, we identify the statistical averages 
of the transit times and the mobilities from the transit time 
and mobility distributions, and match them against the transit 
times and mobilities determined from the geometric analysis 
of photocurrent transients. 

Although convenient to use, it is, however, well known 
that the geometric method extracts single value transit times 
which account for the fastest carriers only. 34 ' 35 We demon- 
strate below, that the corresponding mobilities are not related 
to the average values and that the relative deviation depends 
strongly on temperature and electric field. Furthermore, the 
geometric method is more prone to evaluation errors, as, in 
the case of systems exhibiting a Gaussian DOS, the fitted re- 
gions just span over a very small part of the current transient 
on a double-logarithmic scale and the regions' boundaries are 
judged by the eye, which may lead to huge uncertainties. On 
the other hand, the Scott method, as will be shown below, 
needs just one fit spanning over more than one order of mag- 
nitude in time. Novikovef al. explicitly discussed the conse- 
quences of using poorly defined transit times or mobilities in 
Ref. 36 and asked for a safe and reliable procedure for the 
analysis of ToF data. We believe that the new approach pre- 
sented here fulfils this request, by experimentally providing 
detailed and comprehensive information, as it is usually only 
known from Monte Carlo simulations. 



II. EXPERIMENTAL METHOD 

The studied devices were prepared in a diode configura- 
tion, with P3HT 4002E (~ 94% regioregular; Rieke Metals 
Inc., without further purification) as the active material, drop 
cast on ITO covered glass with an aluminium electrode ther- 
mally evaporated on top. The sample thickness and excitation 
intensity were chosen in a way to ensure a photogeneration 
of charge carriers confined to the first 10% of the bulk thick- 
ness. For the optical excitation we used the second harmonic 
of a pulsed Nd:YAG laser emitting at X = 532nm, close to the 
absorption maximum of the polymer, 37 together with optical 
intensity attenuators. 

Adjusting the polarity of the applied constant electric field, 
one can select the type of the charge carriers being dragged 
through the bulk and extracted at the counter-electrode. Due 
to the extraction, the photocurrent steeply drops towards zero 
until no mobile charge carriers are left within the device. 
However, the photocurrent decreases even before extraction 
due to energetic relaxation of the charge carriers within the 
DOS. 16 " 18 ' 21 Associated with this relaxation is a dispersion 
of the charge carrier package being stronger than predicted 



for purely Gaussian (non-dispersive) transport, together with 
a decreasing mean drift velocity with time. 16 ' 38 As a conse- 
quence, the distribution of transit times is broadened and thus 
the drop in the photocurrent at the extraction is blurred. 

The ToF transients for holes in P3HT where measured in 
the temperatures range between T = 1 10 — 300 K and at dif- 
ferent electric fields between F = 1.3 • 10 7 - 1.9 • 10 8 L 



III. EVALUATION METHOD 

Before describing the Scott method, it is worth having a 
closer look at the photocurrent decay prior to extraction of 
charge carriers, as its functional approximation is an essential 
part of this method. 

In case of a purely exponential DOS this transient decay can 
be explained within the framework of Scher and Montroll 15 
yielding a power law for the pre-transit (jo(f) « t 
photocurrent. Note that this behaviour in combination with 
the predicted post-transit current j{t) °c f-(!+°w) builds the 
basis of the geometric analysis. 33 In contrast to the exponen- 
tial DOS, where transient measurements strictly follow these 
power laws at least "as long as one would be able to mea- 
sure", 39 a Gaussian DOS yields a saturation of the transient 
photocurrent in time, thus it is a power law with a time de- 
pendent exponent y'o(f) x t a ('\ The time dependence of oc(f) 
is related to experimental parameters, such as the temperature 
and the applied electric field, and material properties (ener- 
getic and spatial disorder). 16 ' 38 Unfortunately, a plain func- 
tional description of a(f) and, thus, the photocurrent decay 
jo(t), similar to the power law relation for a exponential DOS, 
is not available. However, taking into account that a(f) is a 
gradually decreasing function and that lima(f) = 0, we sug- 

gest the following empirical approximation of jo(t ): 

jo(t)=f r a V (la) 
a(f) = ( 7 ) . (lb) 

7q is a scaling factor determining the photocurrent at t = 1 s, 
while x is a lifetime and P is a stretching parameter defining 
the dynamics of the decay. We confirm the accuracy of this 
expression in Sec. V using one-dimensional Monte Carlo sim- 
ulations and show that all three empirically defined fitting pa- 
rameters exhibit a strict temperature dependence. The fidelity 
of Eqn. (1) can already be seen from Fig. 1. 

In order to describe the full transient including extraction, 
the integrated number of charge carriers arriving between t 
and df, p tr {t), is subtracted from jo(t): 

j{t) = k{t)(l- f o Ptr{i)^j . (2) 

Hence, the transit time distribution p tr {t) can be calculated 
from the experimental photocurrent j(t) and jo(t) as 

n M d M W 
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TABLE I. Overview over the four different mobility definitions used throughout this Paper. 



Name relation to transit time 



physical context 



Scheme 



Mtr.m 



Scott method 

1 \ ensemble average mobility; determines 

t tr I the macroscopic photocurrent 

1 related to average transit time; often 

k tr \ used in Monte Carlo simulations 



logG) A 



<t, 1 > <V> 




log(t) 



Mgeo 



Ml/2 



geometric method 



logG) A 



1 



l tr,geo 

i 

t tr.\/2 



time at which the first charge carriers 
reach the counterelectrode 

time at which ~ 50% of the charge car- 
riers were extracted 




I L 



log(t) 



Six exemplary transit time distributions calculated from 
measured and simulated ToF transients are plotted in the 
lower part of Fig. 1. A transit time distribution can 
be transformed into a mobility distribution by: pJji) = 
(F-t 2 /d)p tr (t). 

The comparison of the Scott and the geometric method is 
based on the comparison of four explicit mobility definitions. 
To facilitate the later discussion we shortly summarise the dif- 
ferences between the compared mobilities. In Tab. I we depict 
and list how they are obtained, we assign them to their corre- 
sponding method and put a note onto their physical relevance. 
In general the relation between mobility and transit time is: 
p = hj (t tr ■ F), with the sample thickness L and the applied 
electric field F. Thus, the mobility is commonly determined 
in ToF transients from a particular transit time t tr . The mo- 
bility labels and the relation to the corresponding transit times 
can be found in Tab. I and the transit times are obtained as 
follows. 

In the geometrical method t trgeo is identified from the 
intersection of two power laws fitted to the plateau (j pr e) 
and the trailing edge (j pos t) of the photocurrent transients 
in double-logarithmic presentation, and t tr l / 2 is the time at 
which 1 /2 ■ j pre — j post (see lower scheme in Tab. I). From the 
extracted distributions the mathematical expectations can be 
calculated using standard probability theory, 40 ' 41 which gives 
the relation between the transit time (t tr ) and the transit time 
distribution pt r (t) as 

{ttr) = J t tr p, r (t)dt . (4) 

Consequently, the stochastic al average of the mobility (p), 
called p m is 

Mm = (p)= J PP^ipW 



Note that generally the mathematical expectation of the mo- 
bility (p) is not identical to p trjn , which is determined by 
(t lr ), as the former is proportional to the average of the in- 
verse transit time , while the latter is determined by the 
average transit time t tr . Consequently, p m and p trm may yield 
very different values, as yet the mean mobility p m is mathe- 
matically less affected by the fraction of slow charge carriers 
than ptr.m- The latter mobility is often found in literature as an 
output parameter treating the charge transport in disordered 
materials using Monte Carlo 21 or Master Equation 1 ' 24 simu- 
lations. p m , instead, is linked to the average charge carrier 
transport and determines the macroscopic current density by 
j = e-njp- Pp(p)dp F = enp m F. 

IV. SIMULATION METHOD 

In order to verify our approach, a one-dimensional Monte 
Carlo simulation 42 ' 43 was implemented to study charge carrier 
transport by hopping. It considers the capture and the emis- 
sion of charge carriers according to the multiple trapping and 
release (MTR) model. Before further specifying the simula- 
tion procedure, let us first recall some details of the MTR- 
model developed by Schmidlin etal. 44 and Noolandi etal. 45 
Originally the presence of extended states in the valence and 
conduction bands and localised trapping states in the band 
gap were considered. These states are separated by the mo- 
bility edge and only charge carriers which are thermally ac- 
tivated to extended states above the mobility edge contribute 
to charge transport. In disordered organic semiconductors in- 
stead, the valence and conduction bands are described in terms 
of a Gaussian distribution of localised states. Transport in 
these localised states occurs via hopping. In such a system, 
the transport energy 46 is a more relevant concept to describe 
the hopping transport in terms of the MTR. The transport en- 
ergy E trans is a particular temperature-dependent energy level 
within the DOS. 17,47 The transport resembles the MTR pro- 
cess near and below this energy, by a quasi-free transport until 
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trapping of charge carriers above E tram and a thermal emission 
back to E trans of charge carriers below. 

Hence, in our simulation the release, i.e., the emission from 
a trap state, transfers a charge carrier to the transport energy 
Etrans by thermal activation, and above the transport energy 
Etrans a charge carrier moves quasi-free until it is trapped 
again or reaches an electrode and is extracted. The time de- 
pendence of the model allows the simulation of transient mea- 
surement techniques such as ToF. 

Initially, all charge carriers are set in vicinity to the illumi- 
nated electrode and reside at the transport energy E trans . Dur- 
ing the time until capture 



Experiment 



MC-Simulation 



to 



(6) 



the charge carrier at the transport energy move quasi-free 
with the mobility po{F, T) a distance Ax = X c ■ fio(F,T) ■ F. 
Where To = [vo • exp (— 2 • y- a x )]~ is the trapping rate with 
the attempt-to-escape frequency Vo, the inverse localisation 
length y and the intermolecular distance a x . £ is a uniformly 
distributed random number between ]0, 1], F is the electric 
field and T the temperature. The charge carrier is then ran- 
domly trapped into one of the trapping levels E trap with a 
probability according to the Gaussian energy distribution of 
these levels. The release time from a trap level is given by 



; To ■ exp 



Ftrans F tr ap &x'F 



■ln(C) 



(7) 



where kg the Boltzmann constant. From the temporal evolu- 
tion of the charge carriers' positions in the sample, a mean 
current can be calculated by 



j{t) = e-n Q -v(t) 



(8) 



where e is the elementary charge, «o the number of initially 
generated carriers and v(f ) the carriers' mean velocity. 

By suitable parameterisation of the model, experimental 
measurements were reconstructed by simulation, offering a 
deeper insight into microscopic charge transport phenomena. 
To decrease the number free fitting parameters we firstly fixed 
some of them to experimentally determined values as well as 
physically reasonable literature values. The energetic disor- 
der parameter was set to be o = 69meV in accordance with 
the extracted temperature dependence of the zero-field mobil- 
ity in terms of the Gaussian disorder model. 21 The attempt-to 
escape frequency was chosen equal to Vq = 1 • 10 
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TABLE II. Simulation parameters. The parameters were determined 



by reconstruction of experimental ToF transients. 


Parameter 


Value 


Unit 


Source 


o 


69.0 


meV 


A,(l/r 2 ) 21 


Vo 


MO 13 


s- 1 


Ref. 48 


Etrans 


0.0 


cV 


Approx. 


Y 


3.91- 10 9 


m- 1 


Fit 


a x 


1.26- 10~ 9 


m 


Fit 




(0.8-2.15) • 10- 8 


m 2 /Vs 


Fit 



100r 



10 - 



' inc! temperature' 

ee- 




6 - 
4 
2 




inc. field 




inc. temperature 



Ptr(t) 

after eq.(3) 
Ptr(t) 

after eq.(7) 
J I L 




10" 



10" 10"' 
time t / s 



10" 



10" 10"' 
time t / s 



FIG. 1. Upper part: Measured ToF transients of hole currents 
in P3HT at T — 260 K for different applied voltages F = 6.7 • 
10 7 Jj, 1.1 • 10 8 1.5 • 10 8 ^ (left). Simulated ToF transients F = 
1 • 10 7 ^ for different temperatures T = 300K, 250K, 200K (right). 
The pre-transit decay of the photocurrent (dotted lines with markers) 
is reproduced by the fit to Eqn. (1) (lines without markers). Lower 
part: Transit-time distribution ptr(t) for the experimental (left) and 
the simulated (right) transients calculated from Eqn. (3). Note that 
the extracted distributions (solid lines, right) match exactly the distri- 
butions (dotted lines, right) calculated by Eqn. (9) from the temporal 
decrease of the amount N of charge carriers inside the device. 



Ref. 48. We assumed a fixed transport energy at E trans — OeV, 
because we found just little impact of a temperature depen- 
dent transport energy. Further input parameters are the in- 
verse localisation length y and the intermolecular distance a x 
and the charge carrier mobility po(F,T). To fully calibrate 
the simulation we simultaneously fitted a set of measurements 
at different temperatures and electric fields using the global 
minimisation algorithm of differential evolution. 49 The used 
parameters are listed in Tab. II, where we explicitly mark a 
parameter as fixed or fitted. 

Note that from the one-dimensional Monte Carlo simula- 
tions the transit time distribution p tr (t) can be calculated di- 
rectly by the monitored gradual decrease of the number of 
charge carriers inside the device N(t). In this case p tr {t) be- 
comes 



Ptrit) = - 



dN(t) 
dt 



(9) 



V. RESULTS 

A. Verifying the initial photocurrent approximation Eqn. (1) 

At first we demonstrate the quality of the empirical expres- 
sion Eqn. (1) of jo(t). To this end, we simulated a set of tem- 
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perature dependent (T = 200 — 300 K) ToF transients using 
our one-dimensional Monte Carlo simulation and compare the 
transit time distributions obtained by the Scott method with 
the temporal decrease of charge carriers inside the device ac- 
cording to Eqn. (9). The parameterised values used in the sim- 
ulation are listed in Tab. II. In the upper right part of Fig. 1 the 
fit of Eqn. (1) to the simulated pre-transit photocurrent de- 
cay is plotted together with the full simulated transients. The 
approximation of the photocurrent by jo(t) is excellent over 
four orders of magnitude in time. The equivalence of the ex- 
tracted distributions is shown in the lower right part of Fig. 1 : 
both transit time distributions match perfectly well over the 
whole studied temperature range and verify the validity of the 
extrapolation of 70(f) to post-transit times. The correspond- 
ing average transit times of the simulated ensemble of charge 
carriers and the values extracted from the transit time distribu- 
tions p tr {t) by Eqn. (4) are equal within the range of numerical 
precision (A (t tr ) / (t tr ) < 0.5%). 

Further evidence for Eqn. (1) comes from the fact, that al- 
though Eqn. (1) is an empirical approximation the three fitting 
parameters follow a strict temperature dependence, as can be 
seen in Fig. 2. Before we rate these dependencies, let us first 
have a deeper look at the mathematical meaning of jL x and 
P and their relation to physical properties. 

T defines the temporal current saturation. Due to the defini- 
tion of Eqn. (1) jo(t ) has two fixed current values independent 
of the parameter p, /'o(f = 1 s) = yjj and /'o(f = x) = j% -x. Con- 
sequently, as long as the current saturation time X is shorter 
than 1 s j'q is should be proportional to the equilibrium mobil- 
ity. According to Nikitenko etal. 16 this should be the case for 
(a/ (kb -T)) 2 < 5.4. Finally, p specifies the curvature between 
these two fixed points in time at t = X and t = 1 s — the larger 
P the faster the decay. Thus, as relaxation is accelerated with 
increasing temperature, p increases with temperature. 

We found that j'q and the lifetime x scale as 



Jo x exp 



X 0= exp 



-0.45 • 



-0.67 ■ 




respectively. The exponent P depends on 



k R T 



p = 0.64- 



k n -T 



0.031. 



Therefore, for the simulated temperature range in Fig. 2, 
(a/ (kb ■ T)) 2 < 4 Jq should be proportional to the equilibrium 
mobility. And, indeed, the prefactor in the exponential func- 
tion of j'q is found to be 0.45, very close to the predicted value 
of 0.5 16 ' 50 in case of a constant transport energy E trims = 0. 
X and p have very similar prefactors, but inverse temperature 
scaling. However, the physical origin of this scaling might 
be the same, as both parameters define the dynamics of the 
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FIG. 2. Temperature dependence of the empirical fitting parameters 
jn, x and p together with exponential and linear fits. log(jQ) and 
log(t) are proportional to l/T 2 , while (3 scales linearly with T 2 . 



relaxation. Thus, it is very likely that they are related to the 
average charge carrier trapping as well as the charge carrier 
release time limiting relaxation. 2,16 

Based on these results, we conclude that Eqn. (1) is a suit- 
able approximation for the photocurrent decay and that, thus, 
the Scott method can be applied to systems with a Gaussian 
DOS. 



B. Evaluating hole photocurrents with the new approach 

We examined ToF transients measured on P3HT utilis- 
ing the presented approach. We compare the outcome with 
the results obtained by the conventional geometric analysis in 
Sec. VD. In the upper left part of Fig. 1 three ToF-transients, 
parametric in the applied electric field at T = 260 K, are plot- 
ted together with fits of 70(f) according to Eqn. (1) to the 
initial decay of the experimental photocurrent j(t). The ex- 
cellent agreement between the data and the 70(f) over almost 
two orders of magnitude in time was achieved for all mea- 
sured temperatures and electric fields. However, it is ex- 
pected that at very short times 70(f) is not able to approxi- 
mate the photocurrent j(t), since then the current decay on 
the one hand involves a contribution of the quickly extracted 
electrons at the illuminated contact and on the other hand 
is limited by the setup's RC time. The respective transit 
time distributions, plotted in the lower left part of Fig. 1, 
were calculated from the experimental current density j(t) as 
well as the fitted decay 70(f) according to Eqn. (3). These 
transit time distributions — and consequently the mobil- 
ity distributions — are approximately inverse Gaussian 51 " 53 
distributed, in good agreement with distributions obtained 
from measurements on poly(2-methoxy,5-(2'-ethyl-hexoxy)- 
p-phenylene vinylene) (MEH-PPV) thin films performed by 
Rappaport etal. 20 ' 31 and on diketopyrrolopyrrole-naphthalene 
copolymer (PDPP-TNT) thin film transistors by Ha etal. 54 In- 
teresting to note is that inverse Gaussian distributions are well 
known to describe the "First Passage Time Probability" 52 ' 53 
since Schrodinger's theoretical work about brownian motion 
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FIG. 3. Upper part: Mobility distributions p^ijx) for three different 
temperatures T = 130, 190, 300K (from left to right) at a fixed elec- 
tric field of. F = 1.2- 10 s jj. The vertical straight lines illustrate the 
four different definitions of mobility fi trm (dash-double-dotted), jj m 
(dashed), fi geo (solid) and jjijo (dash-dotted) for T = 130K (left) and 
T = 300K (right). Lower part: Conductivity distributions p iU (^) • fJ 
corresponding to the distributions above. 



of particles under a constant force causing drift. This proba- 
bility is actually equivalent to the distribution of transit times 
considered here if a non-dispersive Gaussian charge carrier 
transport, which is nothing other than a brown motion with 
drift, is assumed. 

Derived mobility distributions p^(/j) are displayed in the 
upper part of Fig. 3 for three different temperatures at a fixed 
electric field. Note the different asymmetric shape of p^{p) at 
different temperatures: For the high temperature (T — 260 K), 
the maximum of the distribution is shifted to higher mobili- 
ties, according to many fast charge carriers, whereas the left 
edge is slowly decreasing towards slower mobilities. For the 
intermediate temperature (T = 190 K), the shape is quite sym- 
metric, but broader, with an additional long tail of slow charge 
carriers. At low temperature (T = 140 K) the distribution's 
maximum shifts to lower mobilities and the width is further 
broadened. 

A more intuitive representation of the extracted distribu- 
tions is given by the conductivity distribution /^(/j) • /j plotted 
in the lower part of Fig. 3. It shows the contribution of each 
mobility to the total current density j t . From Fig. 3 it can 
be seen that, although an enhanced number of low mobility 
charge carriers show up in the left edge of the mobility dis- 
tribution at T = 140K, their contribution to the total current 
density j t is small as seen from the conductivity distribution in 
the graph Fig. 3 below. Actually the main contribution comes 
from charge carriers with a mobility around fi m . We point out 
that, using the transport energy definition from Schmechel, 23 
H m is the mobility of charge carriers around the transport en- 
ergy. 



FIG. 4. Upper part, Left axis: Fraction of charge carriers, which 
are extracted before the geometric transit time t lr ^ eo , to all charge 
carriers extracted. A temperature and field dependent deviation is 
seen by the straight lines, which were added as a guide to the eye. 
Upper part, Right axis: Fraction of charge carriers, which are ex- 
tracted before /„-i/2, to all charge carriers extracted. Lower part: 
Temperature and field dependence of differently defined charge car- 
rier mobilities. As the geometric transit time only accounts for the 
fastest charge carriers, the mobility at the geometric transit time fj geo 
lies at systematically higher values. The mobility fj trM from the mean 
transit time together with its standard deviation is indicated by the er- 
ror bars. The mean mobility jA m is physically most relevant for the 
description of an experiment. All mobilities show similar field and 
temperature behaviour. 



C. Comparison of the new approach with the common 
geometrical method 

Next we reveal the advantage and accuracy of the Scott 
method by comparing the results of our data evaluation 
on experimental data with results obtained with the well- 
established geometric approach. Therefore, we remind the 
reader of the two different mobility definitions for each of the 
two approaches listed in Tab. I. 

The four different mobilities are compared in Fig. 3 for two 
temperatures. For both temperatures fi geo has the highest value 
followed by /j m and /j trm is further shifted to lower mobilities. 
/j 1/2 is always situated between /j m and /j„- m . For the two cases 
shown in Fig. 3 the relative error of the geometric mobility 
to the mean mobility is r| = 0.10 at T = 300 K; it strongly 
increases at T = 130K to r| = 0.64. 

In order to quantify the deviation of the geometric transit 
time from the physically more relevant transit time distribu- 
tion, the normalised density of charge carriers extracted until 
the transit time t trgeu is studied in the upper part of Fig. 4. 
While the geometric method accounts for the fastest 40% of 
charge carriers at room temperature and high fields, this value 
drops to below 5% at lower temperatures and fields. Long 
since it is known that p. geu accounts only for the fastest charge 
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carriers 34 ' 35 and, thus, it was often recommended to use /Ui/2 
instead. 35 However, the normalised density of charge carri- 
ers extracted until shows also a slight field dependence 
at room temperature varying between 50% — 60%. A strong 
negative field dependence is instead detected for lower tem- 
peratures. However, at high electric fields nearly no tempera- 
ture dependence is obtained and ~ 60% of the charge carriers 
were extracted until t tr ^i%. 

As the actual charge transport is determined by the mean 
mobility, a systematic field and temperature dependent error 
is made when relying on the geometric method. Although, 
the mobility values determined from t tr ^n better represent the 
mean values, for very dispersive transients with a bad signal- 
to-noise ratio f rr i/2 can not be obtained from experimental 
data, while a transit time distribution can still be extracted. 



D. Extracting an accurate energetic disorder parameter o for 
hole transport in P3HT 

A further advantage of the Scott method is that besides the 
mathematical expectations the mobility and transit time can 
be provided with their standard deviation as a measure of the 
transport dispersion. In the lower part of Fig. 4 we compare 
the field dependence of the different mobilities pt geo , fi m and 
H trm (pin is left out for clarity). Additionally we plotted the 
standard deviations of the mean transit time t trm accounting 
for the transit times shorter than t trm and for the longer tran- 
sit times, respectively. As in Fig. 3, Fig. 4 confirms that the 
geometric mobility /j^,,, always exceeds both /j m and fJ trM . 
Furthermore, the geometric mobility fi ge0 m hardly lies within 
the bounds of the mean transit time mobilities' standard devi- 
ation. Only at high electric fields and high temperatures this 
is the case. It can be nicely seen that the standard deviation 
of the short transit times in Fig. 4 shows a negligible field, but 
strong temperature dependence, as in Ref. 21 Yet, the elec- 
tric field and the temperature strongly increase the standard 
deviation for the long transit times. We attribute this strong 
asymmetry of the transit time distribution ptr(t) around t tr ,mto 
field assisted diffusion. 2 ' 55 ' 56 A detailed investigation of this 
phenomena will be published elsewhere. 

Finally, in terms of the framework of the GDM, 21 we can 
use of the extracted temperature and field dependent ToF mo- 
bilities to gain the characteristic energetic disorder a of the 
DOS. Bassler etal. used the mobility given by the mean tran- 
sit time iutr,m for the derivation of their famous parametric 
GDM mobility equation [Eq. (15) in Ref. 21]. Consequently, 
as our new approach allows for the extraction of this quantity, 
we can calculate an accurate a for the hole transport in P3HT. 
To this end, we demonstrate in Fig. 5 the different scaling 
of the zero field geometric mobility fJo,geo and the zero field 
mean transit time mobility Hoiftr,m) with temperature. The 
extracted energetic disorder O geo t, — 59meV is significantly 
lower by about 14% for the geometric mobilities compared 
to O m j, = 69meV calculated from pL t r,m- The found system- 
atic error of the geometric method leads, thus, to a wrong 
estimate of the energetic disorder parameter a. Comparable 
values for G geo h were found by Mauer etal. 51 using the geo- 
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FIG. 5. Zero field mobilities for geometric /J0,geo an d mean tran- 
sit time mobilities /JQj r . m ). The temperature dependent error of the 
geometric mobility can be clearly seen and results in an underestima- 
tion of the extracted energetic disorder oof the DOS. The transition 
to a weaker temperature dependence at low temperatures is more at- 
tributed to a violation of the common Boltzmann approximation than 
to the transition from dispersive to non-dispersive transport. 



metrical method. For P3HT Sepiolid P100 and Sepiolid P200 
(BASF SE) o = 56meV and o = 58meV was obtained, re- 
spectively. However, our results suggest that these values are 
also underestimated and we a propose a « 64 me V for the hole 
transport in Sepiolid P100 and a « 66meV for Sepiolid P200. 

Besides the impact of the temperature and field dependent 
error of fJo,geo, its deviation to the mean mobility /j m can 
also be understood regarding the energetic distribution of the 
charge carriers. While the fastest charge carriers just hop via 
states around the transport energy and thus "feel" a smaller en- 
ergetic disorder, the mobility calculated from the distribution 
of all charge carriers accounts for all states in the DOS. 

Note that, besides the extracted statistical parameters dis- 
cussed herein, the mobility and transit time distributions in- 
volve much more information about the energy distribution of 
the DOS. At temperatures below T w 190K the slope of the 
zero field mobilities in Fig. 5 changes and become smaller. 
The slope of the geometric values is also lower and the rela- 
tive difference is 11%. Albeit the evaluated values of C geo j = 
3 1 meV and a,„j = 35 meV are not related to the energetic 
disorder a of the DOS and are just given to number the dif- 
ferent slopes. Furthermore, we do not necessarily relate this 
to a transition from dispersive to non-dispersive transport, as 
commonly suggested. First of all, we do not see a significant 
change in the photocurrent shape around T w 190K and sec- 
ondly we observe the very same temperature behaviour of the 
mobility from steady-state field effect measurements 58 on the 
same material system. As, in the case of the Scott method, we 
can rule out a systematical error due to the evaluation proce- 
dure as the nature of this behaviour, we attribute this more to a 
violation of the common Boltzmann approximation utilised in 
the GDM, which predicts the constant slope of ln(/j) vs 1/T 2 . 
Instead, at low temperatures the full Fermi-statistics need to 
be taken into account considering the charge carrier density 
effects on the mobility. 23 ' 25 
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VI. CONCLUSION 

We have introduced a new approach to derive the charge 
carrier mobility and transit time distribution in energetically 
disordered semiconductors exhibiting a Gaussian DOS from 
photocurrent transients. 

It is known that the initial photocurrent decay in these sys- 
tems cannot be described by a plain analytical relation, 16 ' 38 
such as the power law decay in systems with an exponential 
distribution of localised states within the band gap. 15 How- 
ever, a prerequisite of this new approach, which is based on 
the earlier work of Scott etal., 30 is the functional description 
of the photocurrent decay in the pre-transit regime and the ex- 
trapolation into the post-transit regime. Here, we proposed 
an empirical approximation [Eqn. (1)] for the pre-transit pho- 
tocurrent, which accurately fits our experimental data mea- 
sured on P3HT. The validity of this functional approximation 
and the reliability of its extrapolation were shown by one- 
dimensional Monte Carlo simulations. The transit time dis- 
tributions, which we either determined by evaluating simu- 
lated photocurrents with the new method or directly from the 
temporal decrease of the amount of charge carriers inside the 
simulated device, match exactly. Furthermore, the parameters 
of Eqn. (1), although empirical, were found to follow a phys- 
ically reasonable temperature dependence. 

By applying this method to hole ToF transients measured 
on P3HT over a wide temperature and electric field range we 
demonstrated its convenient usage and the superiority com- 
pared to the commonly used geometric method. In contrast to 
the latter the new approach has the following advantages: 

• able to extract parameter distributions instead of a sin- 
gle value only, providing detailed and comprehensive 



information, as it is usually only known from Monte 
Carlo simulations 

• the extracted distributions provide detailed informa- 
tion about the transporting medium, such as, e.g., a 
statistical mean mobility determining the macroscopic 
photocurrent or the dispersion of mobilities, exhibiting 
much potential for further careful studies and input for 
theoretical considerations 

• less error prone, as the pre-transit fit spans over more 
than one order of magnitude in time and the extracted 
stochastical values are numerically averaged 

Furthermore, we found that the geometric method clearly 
overestimates the mobility defined in terms of a mean drift, 
which is the expectation value [i m of the mobility distribution. 
This deviation becomes even more prominent at lower electric 
fields and temperatures. The analysis of hole photocurrents 
with the new approach revealed an accurate energetic disor- 
der parameter O/, = 69 meV for P3HT. This value was calcu- 
lated from the mobility given by the mean transit time Htr,m, 
as this is the relevant parameter considered in the framework 
of the GDM. 21 By using the less accurate geometrical method 
we determined a ~ 14% smaller value of Oh tgeo — 59 meV 
instead. 
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